Search for: All records

Abstract Earthquake focal mechanisms provide crucial information about subsurface fault geometry and stress orientations. Focal mechanisms are typically inferred through analysis of seismic radiation patterns, for example, using P-wave first-motion polarities, potentially in combination with S/P amplitude ratios, to identify nodal planes. The motivation for this procedure is well-founded, as P- and S-wave radiation patterns depend fundamentally on the fault orientation. However, in practice, S/P amplitude ratio measurements can be strongly influenced by factors that are unrelated to the source mechanism. In this study, I characterize several underappreciated issues with S/P amplitude ratio data that are relevant to focal mechanism inversion. The analysis combines synthetic tests with new waveform measurements from ∼64,000 ML≥1.0 earthquakes in Nevada and California. Key findings include that (1) the statistical distribution of S/P amplitude ratio data differs markedly in shape and width from the theoretical expectation, (2) S/P amplitude ratios decay systematically with source-station distance beyond ∼60 km or so, (3) this distance effect is more severe for smaller earthquakes than for larger ones, and (4) modifying the frequency band in which amplitudes are measured can shift the observed amplitude ratio distribution but does not significantly mitigate issues (1)–(3). Taken together, these findings indicate that S/P amplitude ratio measurements are influenced by differential path attenuation and signal-to-noise effects that are not accounted for with existing workflows. Using independent moment tensor solutions, I systematically test various strategies to incorporate S/P amplitude ratios into focal mechanism solutions. The best-performing strategies transform S/P amplitude data to better match the theoretical expectation. Overall, S/P amplitude ratio data appear helpful in improving a typical mechanism solution, but even with the best-performing strategies considered here, the inclusion of S/P amplitude ratio data is expected to hinder rather than improve the solution for a subset of events. 

The size, frequency, and geographic scope of severe wildfires are expanding across the globe, including in the Western United States. Recently burned steeplands have an increased likelihood of debris flows, which pose hazards to downstream communities. The conditions for postfire debris‐flow initiation are commonly expressed as rainfall intensity‐duration thresholds, which can be estimated given sufficient observational history. However, the spread of wildfire across diverse climates poses a challenge for accurate threshold prediction in areas with limited observations. Studies of mass‐movement processes in unburned areas indicate that thresholds vary with local climate, such that higher rainfall rates are required for initiation in climates characterized by frequent intense rainfall. Here, we use three independent methods to test whether initiation of postfire runoff‐generated debris flows across the Western United States varies similarly with climate. Through the compilation of observed thresholds at various fires, analysis of the spatial density of observed debris flows, and quantification of feature importance at different spatial scales, we show that postfire debris‐flow initiation thresholds vary systematically with short‐duration rainfall‐intensity climatology. The predictive power of climatological data sets that are readily available before a fire occurs offers a much‐needed tool for hazard management in regions that are facing increased wildfire activity, have sparse observational history, and/or have limited resources for field‐based hazard assessment. Furthermore, if the observed variation in thresholds reflects long‐term adjustment of the landscape to local climate, rapid shifts in rainfall intensity related to climate change will likely induce spatially variable shifts in postfire debris‐flow likelihood. 

Abstract The Rock Valley fault zone in southern Nevada has a notable history of seismic activity and is the site of a future direct comparison experiment of explosion and earthquake sources. This study aims to gain insight into regional tectonic processes by leveraging recent advances in seismic monitoring capabilities to elucidate the local stress regime. A crucial step in this investigation is the accurate determination of P-wave first-motion polarities, which play a vital role in resolving earthquake focal mechanisms of small earthquakes. We deploy a deep learning-based method for automatic determination of first-motion polarities to vastly expand the polarity dataset beyond what has been reviewed by human analysts. By the integrating P-wave polarities with new measurements of S/P amplitude ratios, we obtain robust focal mechanism estimates for 1306 earthquakes with a local magnitude of 1 and above occurring between 2010 and 2023 in southern Nevada. We then use the focal mechanism catalog to examine the regional stress orientation, confirming an overall trans-tensional stress regime with smaller scale complexities illuminated by individual earthquake sequences. These findings demonstrate how detailed analyses of small earthquakes can provide fundamental information for understanding earthquake processes in the region and inform future experiments at the Nevada National Security Site. 

The state of Nevada is home to one of the most seismically active regions in the world, with crustal deformation associated with the Walker Lane transitioning into Basin and Range tectonics as one traverses from west to east across the state. Despite hosting numerous prominent earthquake sequences over the past century and beyond, at present, there exists no unified research-quality earthquake catalog for the state and its surrounding region. Here, we present a newly compiled, high-precision catalog of more than 180,000 earthquakes occurring around Nevada from 2008 to 2023. The data processing workflow to create this catalog includes an absolute location step that accounts for topography and 3D variations in subsurface wavespeed, and a relative relocation step that refines event positions using differential times measured from waveform cross-correlation. We also provide an update to the local magnitude scale that better accounts for the observed distance attenuation of waveform amplitudes as well as local site effects. We describe some fundamental insights that can be derived from the new catalog, including regional variations in event depth distributions and sequence clustering statistics, and publish the catalog to the wider community to facilitate future research efforts. 

Abstract Quantifying the size of earthquakes is a foundational task in seismology, and over the years several magnitude scales have been developed. Of these, only scales based on seismic moment or potency can properly characterize changes in event size without saturation. Here, we develop empirical potency–magnitude scaling relations for earthquakes in the western United States, allowing us to translate instrumental magnitude estimates into uniform measures of earthquake size. We use synthetic waveforms to validate the observed scaling relations and to provide additional insight into the differences between instrumental and physics-based magnitude scales. Each earthquake in our catalog is assigned a clustering designation distinguishing mainshocks from triggered seismicity, along with a potency-based magnitude estimate that is comparable to moment magnitude and that can be easily converted into other magnitude scales as needed. The developed catalog and associated scaling relations have broad applications for fundamental and applied studies of earthquake processes and hazards. 

Foreshocks are the most obvious signature of the earthquake nucleation stage and could, in principle, forewarn of an impending earthquake. However, foreshocks are only sometimes observed, and we have a limited understanding of the physics that controls their occurrence. In this work, we use high-resolution earthquake catalogs and estimates of source properties to understand the spatiotemporal evolution of a sequence of 11 foreshocks that occurred ~ 6.5 hours before the 2020 Mw 4.8 Mentone earthquake in west Texas.  Elevated pore-pressure and poroelastic stressing from subsurface fluid injection from oil-gas operations is often invoked to explain seismicity in west Texas and the surrounding region. However, here we show that static stresses induced from the initial ML 4.0 foreshock significantly perturbed the local shear stress along the fault and could have triggered the Mentone mainshock. The majority (9/11) of the earthquakes leading up to the Mentone mainshock nucleated in areas where the static shear stresses were increased from the initial ML 4.0 foreshock. The spatiotemporal properties of the 11 earthquakes that preceded the mainshock cannot easily be explained in the context of a preslip or cascade nucleation model. We show that at least 6/11 events are better classified as aftershocks of the initial ML 4.0.  Together, our results suggest that a combination of physical mechanisms contributed to the occurrence of the 11 earthquakes that preceded the mainshock, including static-stressing from earthquake-earthquake interactions, aseismic creep, and stress perturbations induced from fluid injection.  Our work highlights the role of earthquake-earthquake triggering in induced earthquake sequences, and suggests that such triggering could help sustain seismic activity following initial stressing perturbations from fluid injection. 

Abstract Earthquakes are clustered in space and time, with individual sequences composed of events linked by stress transfer and triggering mechanisms. On a global scale, variations in the productivity of earthquake sequences—a normalized measure of the number of triggered events—have been observed and associated with regional variations in tectonic setting. Here, we focus on resolving systematic variations in the productivity of crustal earthquake sequences in California and Nevada—the two most seismically active states in the western United States. We apply a well-tested nearest-neighbor algorithm to automatically extract earthquake sequence statistics from a unified 40 yr compilation of regional earthquake catalogs that is complete to M ∼ 2.5. We then compare earthquake sequence productivity to geophysical parameters that may influence earthquake processes, including heat flow, temperature at seismogenic depth, complexity of quaternary faulting, geodetic strain rates, depth to crystalline basement, and faulting style. We observe coherent spatial variations in sequence productivity, with higher values in the Walker Lane of eastern California and Nevada than along the San Andreas fault system in western California. The results illuminate significant correlations between productivity and heat flow, temperature, and faulting that contribute to the understanding and ability to forecast crustal earthquake sequences in the area. 

Abstract Understanding the connection between seismic activity and the earthquake nucleation process is a fundamental goal in earthquake seismology with important implications for earthquake early warning systems and forecasting. We use high-resolution acoustic emission (AE) waveform measurements from laboratory stick-slip experiments that span a spectrum of slow to fast slip rates to probe spatiotemporal properties of laboratory foreshocks and nucleation processes. We measure waveform similarity and pairwise differential travel-times (DTT) between AEs throughout the seismic cycle. AEs broadcasted prior to slow labquakes have small DTT and high waveform similarity relative to fast labquakes. We show that during slow stick-slip, the fault never fully locks, and waveform similarity and pairwise differential travel times do not evolve throughout the seismic cycle. In contrast, fast laboratory earthquakes are preceded by a rapid increase in waveform similarity late in the seismic cycle and a reduction in differential travel times, indicating that AEs begin to coalesce as the fault slip velocity increases leading up to failure. These observations point to key differences in the nucleation process of slow and fast labquakes and suggest that the spatiotemporal evolution of laboratory foreshocks is linked to fault slip velocity.